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ABSTRACT 

The spectral-energy (and luminosity) correlations in long Gamma Ray Bursts are being hotly 
debated to establish, first of all, their reality against possible selection effects. These are best 
studied in the observer planes, namely the peak energy S'peak ^^ '■^^ fluence F or the peak flux 
P. In a recent paper (Ghirlanda et al. 2008) we started to attack this problem considering all 
bursts with known resdhift and spectral properties. Here we consider instead all bursts with 
known E°^l^, irrespective of redshift, adding to those a sample of 100 faint BATSE bursts 
representative of a larger population of 1000 objects. This allows us to construct a complete, 
fluence limited, sample, tailored to study the selection/instrumental effects we consider. We 
found that the fainter BATSE bursts have smaller E°^l^^ than those of bright events. As a 

consequence, the E°^l^^ of these bursts is correlated with the fluence, though with a slope 
flatter than that defined by bursts with z. Selection effects, which are present, are shown not 
to be responsible for the existence of such a correlation. About 6% of these bursts are surely 
outliers of the i^pcak ^ ^-iso correlation (updated in this paper to include 83 bursts), since 
they are inconsistent with it for any redshift. -Bp^ak correlates also with the peak flux, with a 
slope similar to the i^poak — iiso correlation. In this case there is only one sure outlier. The 
scatter of the E°^l^-P correlation defined by the BATSE bursts of our sample is significantly 
smaller than the B'poak"'^ correlation of the same bursts, while for the bursts with known 
redshift the Speak — -Eiso correlation is tighter than the Speak — Siso one. Once a very large 
number of bursts with Speak ^^^ redshift will be available, we thus expect that the Speak —Lisa 
correlation will be similar to that currently found, whereas it is very likely that the Speak— Siso 
correlation will become flatter and with a larger scatter. 

Key words: Gamma rays; bursts — Radiation mechanisms; non-thermal — X-rays; general 



1 INTRODUCTION 

The correlation between the peak spectral energy -Bpcak and the 
bolometric isotropic energy Eiao emitted during the prompt (the so 
called "Amati" correlation, Amati et al. 2002) may be a key ingre- 
dient for our comprehension of the physics of Gamma Ray Bursts 
(GRBs). The proposed interpretations explain the -Bpcak — Eiao 
correlation as either due to geometric effects (Eichler & Levinson 
2004; Toma, Yamazaky & Nakamura 2005) or to the radiative pro- 
cess responsible for the burst prompt emission (Thompson 2006; 
Thompson, Meszaros & Rees 2007), though there is no unanimous 
consensus. In addition, there is no agreement about the reality of the 
correlation itself. Indeed, Nakar & Piran (2005) and Band & Preece 
(2005) have pointed out the existence of outliers, while Ghirlanda 
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et al. (2005a) and Bosnjak et al. (2008), considering an updated 
Amati relation, found no new outliers besides GRB 980425 and 
GRB 031203. More recently it has been argued that this correlation 
might be the result of selection effects related to the detection of 
GRBs (Butler et al. 2007 - hereafter B07). 

The existence of an -Bpcak — £"130 correlation was predicted 
(Lloyd, Petrosian & Mallozzi 2000) before the finding of Amati 
et al. (2002). The prediction was based on the existence of a sig- 
nificant correlation in the observer frame between the peak en- 
ergy -E°J^|]j and the bolometric fluence F of a sample of BATSE 
bursts. This finding has been recently confirmed by Sakamoto et 
al. (2008a) using a sample of bursts detected by Swift, BATSE and 
Hete-Il. In particular, they note that X-Ray Flashes and X-Ray- 
Rich satisfy and extend this correlation to lower fluences. Amati et 
al. (2002) discovered the Speak ^ E^ao correlation with a sample of 
12 bursts detected by BeppoSAX with spectroscopically measured 
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redshifts. Later updates (e.g. Lamb, Donaghy & Graziani 2005; 
Amati et al. 2006; Ghiiianda et al. 2008 - hereafter G08) confirmed 
it witli larger samples. In the most recent updates (76 bursts in G08 
and 83 in this paper) the GRBs with measured redshift and -E°c|ij 
define a correlation iJpcak oc Ef^^, with only two outliers (GRB 
980425 and GRB 031203, but see Ghisellini et al. 2006). This last 
sample contains bursts detected by different instruments/satellites, 
i.e. BATSE, BeppoSAX, Hete-II, Konus-Wind (all operative in the 
so-called pre-Swift era) and, since 2005 (mostly) by Swift. These 
instruments have different detection capabilities and also different 
operative energy ranges. 

For these reasons is crucial to answer to the following ques- 
tion: is the -Epcak — ^-iso correlation real or is it an artifact of some 
selection/instrumental effect? 

To investigate this issue we should move to the observer frame 
plane corresponding to the Amati relation (i.e. the E°^^^^-F plane) 
where the instrumental selection effects act. There, we can place 
the "cuts" corresponding to the instrumental selection effects and 
see if the distribution of the data points are affected by these cuts, 
or if, instead, they prefer to lie in specific regions of the plane, away 
from these cuts. 

There are two main selection effects: first, a burst must have a 
minimum flux to be triggered by a given instrument. This minimum 
flux can be converted (albeit approximately) to a fluence by adopt- 
ing an average flux to fluence conversion ratio, as done in G08. 
This is the minimum fluence a burst should have to be detected. 
We call it "trigger threshold", TT for short. Secondly, we need a 
minimum fluence also to find -Bpelk ^rid the spectral shape. In fact 
we can have bursts that, although detectable, have too few photons 
around -Bpel^ to reliably determine E°]^ly. itself. Consider also that 
the limited energy range of any detector inhibits the determination 
of -Bpelk outside that range. We call it "spectral threshold" - ST. 



While both selection effects are functions of E'. 



peak 



we found 



that the latter is dominant for all the detectors. 

In G08 we considered a sample of 76 GRBs with redshifts z 
(which we will call, for simplicity, zGRB sample). These bursts de- 
fine an Amati correlation in the form -Bpeak oc E-^^^ , consis- 
tent with what found by previous works. Moreover, G08 found that 
these bursts define a strong correlation also in the observer frame 

/cobs ^ n,0.39±0.05N 

(-fipcakO^ -r ;. 

We demonstrated that the ST truncation effect is biasing the 
2;GRB sample of Swi/j-bursts (i.e. bursts for which -Ep^ak has been 
determined from the fit of BAT spectra) while this is not the case 
for the no-Swift 2:GRBs (i.e. bursts for which E°]!.ly. has been 
determined from other instruments, namely Konus-Wind, BATSE, 
BeppoSAX and Hete-II) . This leaves open two possibilities: (i) 
if those described above are all the possible conceivable selection 
effects, the no-Swift zGRB sample represents an unbiased sample 
and, therefore, the E°]^lj^-F correlation defined with these bursts 
is real; (ii) there are other selection effects biasing the sample of 
zGRBs. For instance, the optical afterglow luminosity might be 
proportional to the burst fluence, resulting in a bias in favour of 
7-ray bright bursts. Another effect concerns the BATSE bursts, 
that had to be localized by the Wide Field Camera (WFC) of 
BeppoSAX. Although, formally, the TT for the WFC should not 
introduce any relevant truncation, in G08 we have shown that all 
bursts detected by the WFC (with and without redshift or measur- 
able -Bpelk) had fluences much larger than its TT curve. 

To proceed further, in this paper we consider all bursts with 
measured -Bpel^, irrespective of having or not a measured redshift. 

With this enlarged sample we study if the distribution of GRBs 
in the E°^l^-F plane is: (i) consistent or not with the correlation addition to z, also E°°ly. is needed. C07, through the analysis of the 



defined by the zGRBs; and (ii) if it is strongly biased or not by the 
considered selection effects. 

Besides using existing samples of bursts (from Hete-II, 
Sakamoto et al. 2005, from Swift, B07, and from BATSE, Kaneko et 
al. 2006, K06), we collected a Konus-Wind sample from the GCN 
circulars (Golenetskii et al., 2005, 2006, 2007, 2008) and we anal- 
ysed a new sample of BATSE bursts reaching a fluence of 10~^ [erg 
cm" ^ ] . This is the BATSE limiting fluence (ST) which allows to de- 
rive a reliable -Bpjjjj^ from the spectral analysis. We therefore have 
a complete, fluence limited, sample of BATSE bursts. With this we 
can study if there is any E'^'^^^-F correlation and if it is a result of 
selection effects or not. 

By the same token, we can study the correlation between 
Speak and the peak luminosity Liso- This correlation was first 
found by Yonetoku et al. (2004), with a small number (16) of bursts, 
and was slightly tighter than the Speak — Eiso correlation, and had 
the same slope. It will be interesting to see if this is still the case 
considering our sample of 83 zGRBs. We can then investigate the 
instrumental selection effects acting on this correlation by studying 
it in the E'^^-^-P plane, where F is the peak flux. 



2 THE S°^|k-FLUENCE PLANE 

The observer frame E'^^-^-F correlation, found using zGRBs 
(filled symbols in Fig. [TJ, divides the plane in two regions corre- 
sponding to large fluence and low/moderate -Bpjjjj^ (region-I) and 
low fluence and large/moderate -Bpol^ (region-II). The absence of 
bursts in region-I suggests that they are extremely rare because, 
otherwise, they would have been easily detected by present and 
past GRB detectors. The absence of bursts in region-II could be 
due to selection effects. 

In the following we will refer to GRBs with known z as 
zGRBs and we will indicate bursts without measured z as GRBs. 



2.1 GRB samples with redshifts: ^GRBs 

Different detectors/satellites (BATSE-CGRO, Hete-II, Konus- 
Wind, BeppoSAX, Swift) have been contributing the sample of 
GRBs with measured z and -Bp^Jk' possibly introducing different 
instrumental selection effects. G08, by considering the trigger ef- 
ficiency of these satellites (from Band 2003), excluded that this is 
affecting the sample of 76 bursts defining the E°^^j^-F correlation. 
However, a stronger selection effect was studied in G08: to add a 
point in the Speak — Eiso correlation, in addition to detect it, we 
need to determine (through the spectral analysis) the peak energy 
Speak- This requires a minimum fluence. G08 simulated several 
spectra of GRBs by assuming they are described by a Band func- 
tion. The values assumed for the low and high energy indeces are 
fixed to the typical values a = — 1 and /3 = —2.3. By varying the 
peak energy, the fluence and the duration G08 derived the "spectral 
analysis thresholds" ST (shaded curves in Fig. \Q in the E°^l^-F 
plane for BATSE, BeppoSAX and Swift. Details of the simulations 
are given in G08. These curves show that the limiting fluence F 
is a strong function of -Ep^ak- ^ burst on the right side of these 
curves has enough fluence to constrain its peak spectral energy. As 
discussed in G08, Fig.[T]shows that the 27 Swift bursts (filled stars) 
with z and well constrained -Bpoaj^ (from C07) are affected by this 
selection effect (dark and light grey areas labeled Swift in Fig.[Tll. 
Note that this is a small sub-sample of the bursts detected by Swift. 
Indeed, in order to add a point to the Epeak — Eiso correlation, in 
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Figure 1. Distribution in the E°^^y-F plane of tlie GRBs with measured redshift (filled symbols) and bursts without measured z published in the literature 
(open symbols). The bolometric fluence is obtained by integrating the spectrum in the range 1 keV-10 MeV. The bright BATSE sample (from K06) is shown 
by open circles (for well constrained £^°j!|]j) and up/down arrows (when only an upper/lower limit can be set on i?°J^|]j from the spectral analysis of K06). 
Hete-U (Sakamoto et al. 2005) and Swift bursts (B07 - but see text) without redshifts are shown with open triangles and stars, respectively. Shaded regions 
represent the ST curves of minimum fluence, for different instruments (see G08 for more details), down to which it is possible to fit the spectrum and constrain 
the spectral parameters. 



HP^-Swift spectra of bursts with known z, showed that, given the 
limited energy range of this instrument (15-150 keV), the peak en- 
ergy could be determined only for a small fraction of bursts. There- 
fore, these 27 Swift bursts are all the GRBs (up to April 2008) for 
which both the z is known and E'^l^ could be determined from 
the fit of the BAT-Swi// spectrum. The fact that the Swift sample, 
for which the spectral analysis of CO? yielded E'^l^, extends to 
the estimated limiting ST is an independent confirmation of the re- 
liability of the method for simulating these curves. Note that a few 
Swift bursts are below these lines, but in these cases the peak energy 
was found using the combined XRT-BAT spectrum (see G08). 

Instead, pre-5wi/? zGRBs (partly detected by BATSE and 
BeppoSAX - filled circles in Fig. [TJ are not affected by the cor- 
responding ST curves. Note that only for BATSE we could derive 
the ST through our simulations. To this aim, the detector response 
and background model is needed. For Konus-Wind, BeppoSAX and 
Hete-II these informations are not public. However, for BeppoSAX 
we can rescale the BATSE thresholds (see G08 for details). 

The sample of zGRBs considered in this work contains 83 
objects, 76 from the sample collected in G08, plus 7 bursts re- 



cently detected (up to April 2008). For all these 7 bursts the spectral 
parameters come from fitting the Konus-Wind spectra and are re- 
ported in the GCN circulars (see Tab. |A3t . With this updated sam- 
ple we find an Amati correlation with slope s — 0.48 ± 0.03 
and scatter a — 0.24J. The same sample defines also a corre- 
lation (Kendall's r = 0.51) in the observer plane in the form 



e: 



peak 



ocF" 



A way to investigate if the lack of GRBs in region-II of the 
^poak-^ plane, i.e. between the distribution of bursts with z and 
the ST curves defined in GOB, is real or it is due to a still un- 
explained selection effect is to consider all GRBs with well con- 
strained Ep^ly. but without measured z. 



^ The scatter is found constructing the distribution of the logaiithmic dis- 
tance orthogonal to the best fit correlation line, and fitting this distribution 
with a Gaussian. 
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2.2 GRB samples without z 

We consider Hete-11 bursts (Sakamoto et al. 2005), Swift bursts 
(from B07), the bright BATSE sample (K06) and Konus-Wind 
bursts (Golenetskii et al. 2005, 2006, 2007, 2008, GCN circulars). 
Through these samples we populate the E°]^lj^-F plane. 



2.2.1 Bright BATSE bursts 

We have considered the sample of 350 BATSE GRBs published by 
K06. The selected bursts have either a peak photon flux (in the 50- 
300 keV energy range) larger than 10 [photon s~^ cm~^] or a flu- 
ence (integrated above 25 keV) larger than 2x10^'' [erg cm~^]. 
We excluded the 17 events whose spectrum was accumulated for 
less than 2 sec as most likely representative of the short duration 
burst population. With the remaining GRBs, we constructed a first 
sample selecting all bursts whose time integrated spectrum is fit- 
ted with a curved model (Band, cutoff power-law or smoothly- 
joined power law) providing an estimate of -Bpcak- This sample 
contains 279 GRBs. The remaining bursts form another sample pro- 
viding lower/upper limits on -EpeJk- those GRBs fitted with a Band 
or smoothly-joined power law with an high energy photon index 
greater than -2 provide a lower limit, as well as those fitted with a 
single flat power law. On the contrary, bursts fitted by single steep 
power laws (photon index < — 2) provide an upper limit on -Ep^ak- 

Fig. [T] shows the K06 sample (open circles) in the E°^lj^-F 
plane together with the 83 zGRBs. The BATSE ST curves are also 
shown. By comparing BATSE GRBs with 2:GRBs (filled symbols 
in Fig. [Hi we note that the two samples are consistent for -Bpeak 
values in the 100 keV - 1 MeV range. However, note that in the 
K06 sample there are also a few bursts with considerably smaller 
fluence (but similar -Epclk) with respect to jzGRBs. In other words, 
there is an indication of the existence of bursts that lie between 
the limiting BATSE curves and the E°]^lj^~F correlation defined 
by zGRBs (region-II). From Fig. [T] it is clear that the sample of 
bright BATSE bursts is not strongly affected by the corresponding 
ST. However, note that this sample is not appropriate to study this 
issue because it is representative only of very bright BATSE bursts 
and it is not complete in fluence. 

The K06 sample shows a weak -Bpeak"^ correlation with a 
Kendall's correlation coefficient r = 0.13 (3a significance). 



2.2.2 Swift and Hete-II bursts 

Other two samples of bursts with published spectral parameters are 
that of Hete-II and Swift. The two references for the Swift bursts 
are C07 and B07: the former focused on Swift bursts with z and the 
latter considered also bursts without redshifts. The C07 Swi/f bursts 
were included in the sample of the 83 2:GRBs. For the Swift bursts 
without z we consider the analysis performed by B07 (but see also 
Sakamoto et al. 2008b). They analysed GRB spectra with either the 
frequentist method and through a Bayesian method. While the first 
method allows to constrain the spectral E°]^l^^ only if it lies in the 
energy range where the spectral data are (15-150 keV for BAT- 
Swift), the bayesian method infers the peak energy by assuming an 
^poak distribution as prior. For homogeneity with the analysis of 
C07 and with the method used to find the ST, we consider only the 
Swift bursts of B07 without z which have their peak energy esti- 
mated through the frequentist method and for which this estimate 
has a relative error <100%. This choice corresponds to the condi- 
tions of the simulated ST of G08. We found 22 Swift bursts which 
satisfy these requirements. 



The Hete-II group published some spectral catalogs of their 
bursts (Barraud et al. 2003; Atteia et al. 2005; Sakamoto et al. 
2005). Sakamoto et al. (2005) performed the time integrated spec- 
tral analysis of 45 GRBs detected during the first 3 years of the 
Hete-II mission. They performed spectral fits by combining the 
data of the high energy detector (Fregate: 6^00 keV) and the low 
energy coded mask detector (WXM: 2-25 keV). We have consid- 
ered in this sample the 27 bursts without z (the remaining are al- 
ready included in the 2:GRB sample) and whose spectrum is fitted 
by a Band or cutoff power-law model which provides an estimate 

of ^pelk- 

Fig. [T] shows that Swift bursts (open stars) and Hete-II bursts 
(open triangles) are both consistent with the correlation defined in 
the SpJ^lij-F plane by the zGRB sample. Also the Swift sample 
without z confirms the validity of the ST estimates. Note that the 
extension of Hete-II events at very low values of -Bpelk is due to 
the fit of their spectrum with the WXM instrument (see Sakamoto 
et al. 2005). 

Note that in the sample of 83 bursts with z there are also the 
BeppoSAX and Konus-Wind events. No spectral catalog of bursts 
without redshifts has been published to date for these two satellites. 



2.2.3 Konus-Wind bursts 

Preliminary results arising from the fit of Konus-Wind spectra can 
be found in the GCN circulars. We collected a sample of 29 GRBs 
(empty squares in Fig. \Q for which an estimate of -Ep^ak ^fd the 
spectral shape is available and which are not already included in 
the Hete-II or Swift samples considered above. For each burst we 
estimate the bolometric (I-IO'* keV) fluence and the bolometric 
peak flux. The results are listed in Tab. lAll in the Appendix. Since 
Konus-Wind covers an energy range from 20 keV to a few MeV, a 
good spectral analysis of very hard bursts can be performed. Any- 
way, the determination of its TT and ST is not possible, as the back- 
ground model and the detector response are not public. Therefore, 
with respect to the distribution of these GRBs in the E°]^lj^-F plane 
(Fig[T), we can only note that it seems to be very similar to that of 
bright BATSE bursts (empty circles). 



3 FAINT BATSE BURSTS 

Among the above sample, the bright BATSE bursts of K06 is the 
largest and, given the spectral range of BATSE, covers a wide range 
in -EpJ^^j.. However, this sample was selected according to a min- 
imum peak flux or fluence threshold and it is not complete either 
in peak flux and fluence. Furthermore it extends down to a fluence 
larger than the the minimum fluence required to derive E°]^l^^. 

Samples reaching smaller fluences indeed exist: Yonetoku et 
al. (2004) performed the spectral analysis of 745 GRBs from the 
BATSE catalog with flux larger than 2 x 10~^ [erg cm~^ s~^]. 
However, they exclude 56 GRBs for which they find a pseudo red- 
shift greater then 12 or no solution using the -Epcak — iiso rela- 
tion as a distance indicator. Thus the final sample is biased by this 
choice as it most likely excludes the bursts with low fluence and 
intermediate/high peak energy and therefore is not representative 
of the whole sample of BATSE bursts with fluences greater than the 
spectral threshold. 

It is clear from Fig.[T]that for BATSE bursts "there is room" to 
extend the K06 sample to smaller fluences: if the ST for BATSE are 
correct, we should be able to analyze bursts with fluence smaller 
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Figure 2. The E°^^,-F plane for the sample of BATSE bursts. The fluence reported in this plot is the bolometric fluence (1 

are GRBs of the K06 sample with catalog fluence larger than 2xl0~^ erg cm 

correspond to those bursts for which we can only estimate a lower/upper limi 

constrain S°Jj^j. from the spectral analysis. The left and right boundaries of this region are calculated for a burst lasting 5 and 20 seconds respectively (see 

G08 for more details). The dotted line represents the best fit to the combined sample: £^°J^|j. 



correspond to those bursts for which we can only estimate a lower/upper limit to E°^^,. The shaded region represents the minimum fluence requested to 



than 2x 10" [erg cm" ]. Therefore, in order to increase the statis- 
tics and test the density of bursts in region-II we extend the BATSE 
sample to the limiting fluence of 10~® [erg cm~'^]. To this aim we 
selected a representative sample of 100 BATSE bursts with a flu- 
ence F above 25 keV (which is a good proxy for the bolometric 
fluence), within the range 10"^ < F < 2 x 10"^ [erg cm"^]. 
The number of extracted GRBs per fluence bin follows the LogTV- 
hogF distribution and therefore this sub-sample is representative 
of the BATSE burst population in this fluence range (which corre- 
sponds to ~ 1000 events). 



For all these bursts we analysed the BATSE Large Area De- 
tector (LAD) spectral data which consist of several spectra accu- 
mulated in consecutive time bins before, during and after the burst. 
The spectral analysis has been performed with the software SOAR 
v3.0 (Spectroscopic Oriented Analysis Routines), which we imple- 
mented for our purposes. For each burst we analysed the BATSE 
spectrum accumulated over its total duration (which in most cases 
corresponds to the Tgo parameter reported in the BATSE catalog). 



In order to account for the possible time variability of the back- 
ground we modeled it as a function of time (see e.g. K06). 

In most cases we could fit either the Band model (Band et 
al. 1993) or a cutoff power law model. To be consistent with the 
method used to derive the spectral threshold curves of Fig. [T] we 
consider that -Bp^ak i^ reliably determined if its relative error is less 
than 100%. If the relative error is greater or if the best fit model is 
a simple power law we derive the corresponding lower/upper limit. 
In these cases the burst is reported on the plot as an up/down arrow. 

In Tab. IA2I in the appendix we list the bursts of our sample 
together with the results of the spectral fitting. 



3.1 Results for the complete sample of BATSE bursts 

In order to construct a complete sample of BATSE bursts, we cut 
the K06 sample at a fluence F (as reported in the BATSE catalog) 
greater than 2 x 10"^ [erg cm~^] (213 GRBs). This complete sub- 
sample is representative of bright BATSE bursts. To this we add 
the 100 bursts of our representative sample of the 1000 GRBs with 
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(see Fig. |2). The shaded distribution is of the 213 GRBs of K06 and the 
hatched distributions is of the 88 GRBs of our analysis (we have excluded 
upper/lower limits). The fit with gaussian functions of the two distributions 
have scatter cr = 0.18 for the K06 sample and a = 0.29 for our sample. 
The insert shows the scatter of the combined sample (cr = 0.26) (solid line), 
once we take into account that the GRBs analyzed by us are representative 
of 1000 bursts. 



fluences between 10 

T^obs 



and 2xl0~^ [erg cm~^]. Fig.|2]shows the 
pg'Ji. and bolometric fluence (computed in the range 1 keV - 10 
MeV) of the sub-sample of bursts from the K06 sample (open cir- 
cles) together with the 100 bursts of our sample (filled circles). This 
combined sample extends the K06 fluence limit to F > 10~® erg 
CTcT^ . Note that in this figure we plot the bolometric (1-10'' keV) 
fluence estimated accordingly with the best fit model. Its value can 
be different from the fluence value reported in the BATSE catalog. 

The distribution of BATSE GRBs in Fig. |2] defines a corre- 
lation with a large scatter. The Kendall's correlation coefficient is 
r = 0.18 (Jo significant). Since the dimmer part of the burst distri- 
bution in the E'^'^y-F plane is affected by the ST truncation effect, 
we analyzed the correlation also following the method proposed by 
Lloyd et al. (2000). We obtain a Kendall's correlation coefficient 
T = 0.2 (5.5a significant). By fitting with the least square method, 
without weighting for the errors and neglecting the upper/lower 
limits, we obtain Efjl^^ oc _F°^J6±o.o2 (^j^ttg^j ji^g ;„ Fig.|2j. 

In Fig. [3] we show the distribution of the scatter of the two 
samples around the best fit correlation in the Spcak"^ plane. These 
have standard deviation of ct=0. 18 for the K06 sample (solid his- 
togram) and o =0.29 for our representative sample (hatched his- 
togram). The combined sample (solid line in the insert) has a scatter 
distribution with o =0.26. In order to describe the scatter distribu- 
tion of BATSE bursts down to the fluence limit of F > lO"'^ [erg 
cm~'^], we have to consider that our sample of 100 bursts is repre- 
sentative of the entire burst population (a factor 10 larger in num- 
ber) in the fluence range 10"*^ < F <2 x lO"'' [erg cm"^]. We 
fitted the scatter distributions of our (dashed line) and K06 (solid 
line) sample with gaussian functions and combined these best fit 
distributions by renormalizing that of our sample by a factor 10 
(corresponding to the ratio of the extracted bursts with respect to 



Figure 4. The E°\l^ distributions for the K06 GRB sample (213 GRBs - 
shaded histogram) and for the sample that we have analyzed (88 GRBs with 
well determined S°^^j^ over 100 bursts randomly extracted from the BATSE 
LogAf-LogF distribution in the range 10~^ < F < 2 X 10~^ [erg cm"^] 
- hatched histogram). The solid and dashed Unes are the gaussian fits to the 
two distributions. The K-S test (see text) confirms that the shift of E°^^, 

peak 

to lower values for dimmer bursts is statistically significant. In the insert are 
reported these two gaussian fits and their sum (solid line) which has been 
obtained by multiplying the distribution of the 100 GRBs by 10 (dot-dashed 
line) in order to account for the total number of bursts in this fluence range. 



the total number of BATSE bursts in the same fluence bin). The re- 
sult is shown in the insert of Fig.[3](solid line). The combined sam- 
ple (solid line in the insert) has a scatter distribution with a=0.26. 
Fig. |2] shows that there are bursts with low fluence and high 
^poak Slid that the dispersion in -B°j!|ij at low fluence is larger then 
the dispersion at high fluence. However, Fig.|2]also shows that, on 
average, the error on the Fp^ak value increases for smaller fluences. 
This could imply that the larger scatter for smaller fluences is in part 
due to larger errors on F°j!Ji^. A simple way to determine the con- 
tributions to the total observed scatter atot, calculated orthogonally 
to the fitting line, is: 



2 



2 2 

cje cos / 



+ 0-C 



(1) 



where a% is the average relative error on E'^l^, 6 is the angle de- 
fined by the slope of the correlation (whose angular coefficient is 
equal to tan 9 = 0.16) and tJc is the intrinsic scatter of the distribu- 
tion. For fluences greater than 2x10^^ [erg cm~^] (K06 sample) 
o"tot ~ cTc = 0.18 since the errors on FpJ^^j. are small. On the other 
hand, for fluences smaller than 2x10" [erg cm~^] (our sample), 
(je ~ 0.18 and atot ~ 0.29, leading to CTc ~ 0.23, to be compared 
to CTc = 0.18 for fluences greater than 2xl0~^ [erg cm~^]. This 
leads us to conclude that the intrinsic scatter around the best fit line 
increases for smaller fluences. A caveat is in order: the scatter (Tc 
does not take into account lower/upper limits, which also do not 
enter in the derivation of the best fit line. Thus tJc could be larger, 
but only slightly, since the number of upper/lower limits is very 
limited. 

Through our BATSE sample we can also study the Fpjj^j. dis- 
tribution. In Fig.|4]we show the E°^ly, distribution of our BATSE 
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sample (hatched histogram) and that of bright BATSE bursts of K06 
(solid filled histogram - cut at 2x 10~^ [erg cm~^]). 

The shift of the -Ep^ak distribution to lower values for smaller 
fluence selection is statistically significant: the K-S test gives a 
probability P = 5.8 x 10^^^ that the two distributions belong to 
the same population. Similarly to what has been done for the scat- 
ter distribution, we combine the two distributions by accounting for 
the fact that our sample is representative of a larger population of 
bursts in the 10^^- 2x10^^ [erg cm^^] fluence range. The result 
is shown in the insert of Fig. |4] (solid line). This total sample dis- 
tribution has a peak at _B°^ak~160 keV, i.e. smaller than the 260 
keV of bright bursts of K06, and a standard deviation a=0.2S. Al- 
though the widths of the distributions of E°]^ly. can be affected by 
the measurement errors, the central values are not. 



4 OUTLIERS OF THE Speak £iso CORRELATION 

In Fig. [5] we combine, in the E°]^l^^~F plane, our sample of BATSE 
bursts with the zGRBs (solid filled squares) and with the Swift, 
Hete-II and Konus-Wind samples of GRBs without redshifts. 

We note that bursts with known redshift (filled squares) are 
only representative of the large fluence (for any -Epeak) P^t of th^ 
plane. In particular. Fig. [5] shows the existence of bursts with low 
fluence (between F ^ lO"'^ and F ~ 10"^ [erg cm"^]) but E°^ly^ 
larger than 200 keV. These events are not present in the zGRB sam- 
ple. Their absence in the zGRB sample suggests the existence of a 
selection effect. 

However, zGRBs are those defining the Amati correlation (as 
in G08). We do not know if all the other bursts (without redshifts) 
represented in Fig. |5] satisfy this correlation. 

If these GRBs have a similar redshift distribution of those with 
measured z, then it is likely that they would define a rest frame 
Epoak — Eiso correlation with different properties (slope, normal- 
ization and scatter), since some of them stay apart from the -Bpelk" 
F correlation defined by the zGRB sample. On the other hand, 
GRB 980425 and GRB 031203 do have a peak energy and fluence 
consistent with the jzGRB sample, but it is their small redshift to 
make them outliers with respect to the -Bpeak — Eiso correlation. 

The possibility that there is a considerable number of outliers 
of the Epcak — Eiso correlation in the BATSE sample has been dis- 
cussed in the literature (e.g. Nakar & Piran 2005; Band & Preece 
2005 ; K06 - but see also Ghirlanda et al. 2005 ; Bosnjak et al. 2007). 
We can test if a burst is consistent or not with the -Epcak — Eiso 
correlation even if we do not know its redshift. Simply, we as- 
sign to the burst any redshift, checking if there is at least one z 
making it consistent with the correlation. By "consistent" we mean 
the the burst must fall within the 3(7 scatter (assumed gaussian) 
of the correlation. This test was first proposed for the short bursts 
(Ghirlanda et al. 2004a) and then applied to BATSE long GRBs. 
More quantitatively, following Nakar & Piran (2005), we can write 
the Epcak — Eiso correlation as 

f-iTidlF" 



E 



peak 



(1 + 2) = k 



\ 1 + 



Epclk = fcEV(^); fiz) = 



{ATvdlr 



(2) 



(l + z)i+" 

where F is the bolometric fluence. The function f{z) has a max- 
imum (/max) at some redshift and therefore all bursts for which 
Epeak/C^^E") > /max are outliers. We can impose that the con- 
stant k accounts for the scatter of the best fit Epoak — Eiso corre- 
lation, and then find outliers at some pre-assigned number of a. It 



correlation 


sample 


scatter 


K 


s 


-^pcak -'^iso 


zGRB 


0.28 


-18.6 


0.40±0.03 




^GRB (only Swift) 




-11.4 


0.26±0.05 




zGRB (not Swift) 




-20.5 


0.44±0.03 


-E-pcak ~ -Eiso 


zGRB 


0.23 


-22.7 


0.48±0.03 




zGRB (only Swift) 




-16.7 


0.36±0.06 




zGRB (not Swift) 




-24.4 


0.51 ±0.03 


^obs _p 
peak 


zGRB 


0.26 


4.41 


0.39±0.05 




BATSE 


0.20 


3.93 


0.28±0.02 


Eobs p 
peak 


zGRB 


0.23 


4.00 


0.40±0.05 




BATSE 


0.23" 


3.07 


0.16±0.02 



Table 1. Results of the congelation analysis. For each congelation in the rest 
frame and observed plane we give the values of the scatter, normalization 
and slope. The correlations are in the form y = Kx", where y is the loga- 
rithmic observed/rest frame peak energy in units of keV and x is the loga- 
rithm of the luminosity (energy) in erg s~^ (erg) or the peak flux (fluence) 
in erg cm~^ s~^ (erg cm~^). a: this is the value once depurated of the 
contribution to the scatter of the measurement enors (see text). 



is worth to recall that this method assumes that the dispersion of 
points, around the Epeak — Eiso correlation under test, is described 
by a Gaussian function. With this assumption we can state that a 
given GRB is Na inconsistent with the correlation, and quantify 
the probability of having a certain number of outliers lying - say - 
more than 3a away. Since the Amati correlation, as discussed be- 
low, surely incorporates an extra-Poissonian dispersion term (Am- 
ati, 2006), the scatter distribution may not be a Gaussian, but it may 
correspond to the distribution function of this extra term. In other 
words: the scatter of the points around the Amati correlation is not 
due to the errors of our measurements, but reveals the presence of 
an extra-observable not considered in the Amati relation. With this 
caveat, we nevertheless use this assumption (i.e. Gaussianity) for 
simplicity. 

In Fig. |5] the grey area identifies, in the Ep^Jj^-E plane, the 
"region of outliers". Considering only BATSE bursts we can state 
that the 6% of the complete sample considered in this paper is con- 
stituted by bursts which are surely outliers of the Epeak — Eiso 
relation. We can test if these outliers have different spectral proper- 
ties with respect to other bursts (that instead pass the above consis- 
tency test). By comparing their spectral parameters we find that the 
outliers of the Epeak — Eiso correlation have a larger peak energy 
than the total sample of bursts (K-S probability P = 8.7 x 10~^) 
and a slightly harder low energy spectral index a (K-S probability 
P — 10^^). From Fig. [5] we also note that there is no outlier for 
the Epeak-E^ correlation. 



5 THE E, 



peak 



Liso CORRELATION 



Yonetoku et al. (2004, Y04), with a sample of 16 GRBs of known 
z, found that Epoak<x Lfso, where Liso is the isotropic luminosity 
at the peak of the prompt light curve, but calculated using the time 
averaged spectrum (i.e. Epeak and spectral indices), and not the 
spectral properties at the peak flux. 

This correlation appeared to be tighter (but with similar slope) 
than the Epeak — Eiso correlation, as originally found by Amati et 
al. (2002). Since then this correlation has been updated only once 
(Ghirlanda et al., 2005b). 

It is interesting to test if the same conclusions that can be 
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Figure 5. Consistency test of the Spoak — ^iso correlation. The open symbols represent the bursts without redshifts detected by Hete—II (triangles), BATSE 
(open and filled circles), Swift (open stars) as described in the legend. The solid line is the Bpcak — ^iso correlation transformed in the observer E°^^,-F 
plane. The long-dashed line is the 3(t scatter of the Speak — ^iso coiTelation. The "region of outliers" is the grey shaded. Bursts falhng in this region are not 
consistent with the Speak — E\g,a correlation for any redshift: they are outliers at more than Sit (if the scatter distribution is Gaussian, see text). In the upper 
left corned we also show the "region of outliers" of the Speak — S-, correlation (adapted from Ghirlanda et al. 2007) if bursts have a 90° jet opening angle. 
The dotted line is the fit to the zGRB sample (filled squares) and the dashed line is the fit to the complete sample of BATSE bursts described in Sec. 3.1 (see 
also Tab[T) 



drawn for the Speak — E\so correlation (i.e. the presence of selec- 
tion effects and of outliers) can now be extended to the iJpeak — ^iso 
correlation. To this aim we have considered the 2:GRB sample (see 
Tab. |A3| in the Appendix) and we have calculated for all these bursts 
their isotropic equivalent luminosity Liso- This is computed by in- 
tegrating the time averaged spectrum after renormalizing it with the 
peak flux. Note that, strictly speaking, this luminosity does not cor- 
respond to the peak luminosity (see Ghirlanda et al. 2005b), since 
it adopts the time averaged Speak- 

In Tab. lA3l we report the sample of 83 GRBs with their peak 
flux, the energy range where it is computed, the references and Liso- 
To calculate Liso we adopted the same method used to compute 
Siso (see Ghirlanda et al. 2007 for more details). 

In Fig.[6]we show the Speak — Liso and the Speak — Eiso corre- 
lations defined with the 83 GRBs of Tab. I A3 1 The no-Swift bursts 



(empty symbols) and the Swift bursts (filled squares) are shown. 
In both cases, the correlations are highly significant (rank corre- 
lation coefficient are respectively 0.83 with a chance probability 
6.2 X 10"^^ and 0.84 with a chance probability 9.8 x 10~^^). The 
solid lines show the best fit with the least square method (with- 
out accounting for the measurement errors): we obtain Speakoc 
L'^^^o±o.o3 gjjj Ep^^i^tx £;.0-48±003 -p^g gjg jjf jjjg no-5w;/f burst 

sample (dotted line) and of the Swift burst sample (dashed line) are 
also shown. The results of these fits performed considering differ- 
ent samples are shown in Tab.[T] 

Our sample of 83 2:GRBs confirms the finding of Yonetoku 
et al. (2004), even if we obtain a flatter slope. Fitting the scatter 
distribution of the Speak — iiso correlation with a Gaussian we 
derive a = 0.28. Comparing it with the corresponding scatter of 
the Speak — Siso Correlation (a = 0.23) we find that, contrary to 
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Figure 6. -Bpeak ~ ^iso and -Epoak ~ ^iso correlations for 83 GRBs with measured redsliift and spectral parameters. The best fit of the whole sample is 
shown with a solid line. Note that the fit performed on the Swift sample alone (filled squares) has in both cases a veiy fiat slope (dashed line) with respect to 
the slope derived for no-Swiyi" bursts (dotted line). The results of these different analysis are reported in Tab.[T] For an explaination of the flat slope found with 
the Swift sample see G08. 



what initially found by Yonetoku et al. (2004), the scatter of this 
correlation is slightly larger. 

We can investigate if this correlation is affected by any of the 
selections effects that have been studied in G08 for the Bpeak — 
Eiao correlation. In particular we show in Fig.[7]the observer frame 
Spg^jj-P correlation where P is the bolometric peak flux. Note that 
also in this plane the jzGRBs define a strong correlation (dotted line 
- with slope 0.39) and that the GRB samples without z considered 
in this work are consistent with this correlation (differently to what 
happens in the E°]^lj^--F plane). The dashed line represents the best 
fit obtained considering only BATSE bursts. They define a flatter 
correlation (slope 0.28) with respect to the zGRB sample. Note that 
this happens also in the E°]^lj^~F correlation and it is likely due to 
the difficulty of the BATSE instrument to see very low -Bpeak '^t 1°^ 
fluence/peak flux. 

The peak flux P is the quantity on which the trigger condition 
(for most instruments) is determined. We plot in Fig.|7]the trigger 
limiting curves (from G08) as a function of Bpelk- We note that 
for BATSE the TT curve is separated from the distribution of the 
corresponding bursts (open and filled circles). This is because the 
dominant selection effect acting on our BATSE complete sample is 
the ST (see Fig.[2l(. In other words, the bursts that can be displayed 
in the E°^lj^-P plane are not all the bursts that can be detected 
by a given instrument, but only those with a sufficient number of 
photons to make possible the determination of Bp^ak- 

The Hete-Il bursts, instead (triangles), are very near to their 
TT. For this instrument we are not able to determine the ST curves, 
but it is likely that the dominant selection effect acting on Hete-II 
bursts is the need to trigger them. 

For Swift bursts we have an intermediate case: their TT curve 
is not truncating their distribution, even if they lie closer (than 
BATSE bursts) to it. 

Also for the Bpeak ~ Biso correlation we can test if there are 
outliers. Ghirlanda et al. (2005b) tested this through a sample of 



442 GRBs with redshifts derived by the lag-luminosity relation. 
They did not find evidence for outliers. In this work we test the 
Bpeak — iiso Correlation with the same method described above 
for the Bpeak — Biso correlation. In Fig.|7]we show the "region of 
outliers" for the Bpeak — Biso correlation. Only one burst (of the 
K06 sample) is inconsistent with this correlation at more than 3(j. 



6 DISCUSSION AND CONCLUSIONS 

To study the role of possible instrumental selection effects on the 
Amati relation we have focused our attention on the observational 
Bpeak-B plane. Here we can compare the distribution of differ- 
ent samples of GRB (for example, zGRBs and GRB with unknown 
redshift). To this aim we adopt the analysis performed by G08, re- 
ferring to two different instrumental biases: the trigger threshold 
(TT, the minimum fluence derived considering the minimum flux 
required to trigger a burst) and the spectral analysis threshold (ST, 
the minimum fluence needed to constrain the GRB spectral prop- 
erties). These curves depends on Bpeak ^nd define what part of the 
observational plane is accessible. 

First we updated the sample of bursts with redshift, adding 7 
new recent GRBs, for a total of 83 objects. These GRBs define 
a BpeaktxBiso"'** ^'^'^ Correlation in the rest frame, very simi- 
lar to that obtained with previous (and smaller) samples. In the 
observer plane, they define a slightly flatter correlation (BpeakO*^ 
^o.40±o.05-j jjjg scatter of these two correlations is the same (see 
Tab. [Til. As G08 pointed out, the BATSE ST curve is not biasing the 
distribution of BATSE bursts with redshift in the observer plane, 
while the Swift ST could, in the sense that the distribution of Swift 
bursts (with redshift) is truncated by the Swift ST curve. Then why 
the BATSE bursts (with redshift) are not truncated by their corre- 
sponding ST? Is it because of a real, intrinsic correlation or is it 
due to another, hidden, selection effect? One way to answer this 
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Figure 7. Consistency test of the -Epcak — ^iso correlation. Tlie tilled squares are bursts with measured redshifts. The solid line is the Bpoak — ^iso correlation 
transformed in the observer E°^^^-P plane (here it is represented the bolometric peak flux). The long dashed line is the 3a scatter of the Speak — ^iso 
correlation (as discussed in Sec. 5). The shaded triangle delimits the "region of outliers". Bursts falling in this region are not consistent with the Speak — ^iso 
congelation for any redshift: they are outliers at more than 3cr (if the scatter distribution is Gaussian, see text). The dotted and dashed lines show the best fit 
obtained by considering respectively the zGRB sample and the BATSE complete sample. The curves represent the TT estimated for different instruments by 
assuming that the trigger is based on the peak flux criterion. 



crucial question is to analyze all GRBs with -Epeak- even without 
redshift. The BATSE sample of GRBs is the best suited for this aim 
because: i) it contains a large number of bursts; ii) large sample of 
bright GRBs have been already analyzed, and iii) for BATSE we 
already know the ST curve. Then we pushed the spectral analysis 
to the limit, deriving the spectral parameters for a representative 
sample of 100 BATSE GRBs with a (bolometric) fluence between 
10~® [erg cm~^] (corresponding to the ST limit) and 2 x 10^'' 
[erg cm^'^] (the limiting fluence of K06). These 100 GRBs repre- 
sent a large population of 1000 GRBs, in the same fluence range. 
Combining our and the K06 samples we have a homogeneous and 
complete sample, best suited to study how BATSE GRBs populate 
the £^pJj|]j-F plane. Using this complete, fluence limited, sample 
we find: 

• GRBs without redshifts, in this plane, are not spread in the 
region free from instrumental selection effects, but define a corre- 



lation with a flat slope (~ 0.16) and a scatter larger for smaller 
fluences (after accounting for the errors increasing for smaller flu- 
ences). Fig.[8]is a graphic illustration of this: different grey levels 
corresponds to different density of points in the E°^lj^-F plane, 
once we take out the effect of the overall increase in density going 
to smaller fluences (for the LogA'^-LogS' slope). The way we do 
this is the following: we consider different fluence-bins and in each 
of those we count the total number Nf of objects. Then we divide 
this fluence-bin into Bpeak-bins, counting the number of objects 
in each small area, dividing it by Nf . Each small area is then char- 
acterized by a number riij- (between and 1), which corresponds 
to a different level of grey. Note that the data points do not fill the 
entire accessible region of the plane but concentrate along a stripe. 
Note that the shape of this concentration of points is not determined 
by the ST curve, reported in Fig.[8]for a typical burst lasting 20 s. 
The only effect of the ST curve to the found correlation is to cut it 
at the smallest fluences and -Epeak- The very flat slope could be due 
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Figure 8. Graphic illustration of the E°^^, -F coiTelation. We consider 
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different fluence-bins and in each of those we count the total number Nf 
of objects. Then we divide this fluence-bin into E°^^, -bins, countine the 

■* peak ° 

number of objects in each small area, dividing it by Nf . Each small area is 
then characterized by a number rii j (between and 1), which con'esponds 
to a different level of grey. In this way the increasing number of bursts for 
decreasing fluence (the "LogA'^-LogS"' effect) is accounted for, and it does 
not influence n^j. Note that the data points do not fill the entire accessi- 
ble region of the plane (to the right of the ST curve shown in Fig. [2), but 
concentrate along a stripe. 



to the difficulty of having BATSE GRBs with E°^li^ smaller than 
~ 50 keV, whose existence is demonstrated by other instruments. 
However, the paucity of the derived upper limits on iSpclk suggests 
that this effect is marginal. 

• Formally, the scatter is not greater than the scatter of the 
-Epeak — Siso Correlation (both have a = 0.23 once the contribution 
to the scatter of the measurement errors are taken into account). 
Despite that, the entire BATSE sample and the 2:GRB population 
define two E°]^l^^-'F correlations which have significantly different 
slopes. If their redshift distribution is similar, then they will define 
two different correlations also in the Speak — Eiso plane: consider- 
ing then the two samples together, we will define a correlation with 
intermediate slope and a scatter larger than the individual one. 

• If the above point holds (i.e. if the redshift distributions of 
GRBs of unknown redshifts is the same of the zGRBs) then we 
can conclude that there exists an -Bpeak — Siso correlation, not de- 
termined by selection effects, even if its slope and scatter may be 
different from what we know now. We should emphasize that by the 
term "correlation" we mean that GRBs will occupy a "stripe" in the 
-Epeak — £^iso plane with a relatively large scatter (fitting it with x^ 
method one would obtain a very large reduced Xr)- In other words, 
it is very likely that there is another (third) variable responsible for 
the scatter. In fact one finds a tighter correlation considering, as a 
third variable, the jet break time (Ghirlanda, Ghisellini & Lazzati 
2004; Laing & Zhang 2005) or the time of enhanced prompt emis- 
sion (Firmani et al. 2005). Another cause for a large scatter is the 
viewing angle, if a significant number of bursts are seen slightly 
off-axis. 

• In the BATSE sample there are a few bursts with small or inter- 
mediate fluences but large -Ep^ak' not present in the zGRB sample. 
Among them there are some surely outliers of the -Epeak — Eiso 



correlation (as defined by the zGRB sample), i.e. bursts that lie at 
more than 3(t from it, for any redshift. The number of these sure 
outliers is however very small, amounting to the 6 per cent of the 
entire population. 

• We have also investigated the -Epeak — Eiso correlation, and 
its counterpart (E°]!.lj^-P) in the observer plane. First, we partly 
confirm the original findings of Yonetoku et al. (2005, see also the 
update in Ghirlanda et al. 2005b): with the jzGRB sample we find 
a strong Ep^ak — Eiso correlation, whose slope is flatter than orig- 
inally found (s = 0.40 instead of 0.5) and whose scatter is greater 
than the scatter of the Epeak — Eiso correlation. 

• In the observer plane, instead, the Epjj|i.-P correlation of our 
complete sample of BATSE bursts is tighter than the the E°]!.lj^-F 
correlation (a = 0.2 instead of cr = 0.23). Its slope is s = 0.24, 
flatter than the Epeak — Eiso correlation (s = 0.39), but however 
steeper than the Ep^lj^-F slope (s — 0.15). There is only one sure 
outlier. 

• Selection effects are in this case determined by the TT curves. 
These effects are present, being responsible for the cutting at low 
peak fluxes, but they do not influence the slope and scatter for peak 
fluxes larger than the what defined by the TT curves. 

• Considering the jzGRB sample we have that the Epeak — Eiso 
correlation is tighter than the Epeak — Eiso one. Considering our 
complete BATSE sample and moving to the observer plane, we 
have just the opposite: the Epelk-E correlation is tighter than the 



e: 



peak 



-F one. 



• It is then conceivable that the Epeak ~ Eiso correlation, once 
a large number of burst with redshift will be available, will be 
stronger than the Epeak — Eiso one. 

The general conclusion we can draw from our study is that, 
although selection effects are present, they do not determine the 
spectral-energy and spectral-luminosity correlations. These could 
be characterized by a slope and scatter different from what we have 
determined now using heterogeneous bursts samples with mea- 
sured redshift, but we found that Epeak is indeed correlated with 
the burst energetics or peak luminosity. Therefore it is worth to in- 
vestigate the physical reason for this relation. 
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Table Al. Spectral and temporal properties of 29 Konus—Wind bursts without known redshift. The listed fluence and peak flux are estimated in the range 1-10^ 
keV. It is also reported (when available) the estimated burst duration T. The references are: Golenetskii et al., 2005a, 2005b, 2005c, 2005d, 2005e; Golenetskii 
et al., 2006a, 2006b, 2006c, 2006d, 2006e, 2006f, 2006g, 2006h, 20061; Golenetskii et al., 2007a, 2007b, 2007c, 2007d, 2007e, 2007f, 2007g, 2007h, 20071, 
2007J, 2007k; Golenetskii et al., 2008a, 2008b, 2008c, 2008d. The GCN circulars number is reported in the last column. 
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-1.20 




104 ± 19 


1.2±0.3)e-6 


6824 


0.35 ±0.70 


-2.06±0.89 


342 ± 220 


(2.5±2.1)e-6 


3093 


-1.20 




261 ± 225 


6.6±4.3)e-6 


7290 


-0.33±0.25 




1136 ±582 


(3.8±2.0)e-6 


3101 


-1.59±0.22 


-2.01±0.09 


138 ±21 


l.l±0.3)e-5 


7319 


-0.23±0.5 


-2.30 


121 ±81 


(9.6±8.0)e-6 


3177 


-1.09±0.18 


-2.55±0.52 


164 ± 55 


2.1±0.8)e-6 


7374 


-0.81 ±0.42 


-2.26±0.84 


207 ± 142 


(8.6±7.5)e^ 


3217 


-1.60 


-2.20 


231 ±151 


8.0±2.1)e-6 


7387 


0.51 ±1.58 


-2.34±0.21 


83 ±69 


(4.6±4.2)e-6 


3220 


-0.19±0.60 




179 ± 100 


2.7±2.0)e-6 


7504 


-1.34±0.22 


-2.61±0.50 


107 ± 43 


(3.4±1.8)e-6 


3276 


-1.00 


-2.30 


190 ± 53 


5.9±1.8)e-7 


7552 






<30 


(1.5±0.7)e-6 


3319 


-1.02±0.37 




195 ± 134 


l.l±0.9)e-6 


7597 


-1.20 




163 ± 141 


(1.4±1.0)e-6 


3516 


-1.23±0.17 


-2.28±0.66 


313 ±180 


9.4±3.8)e-6 


7638 


-0.66±0.52 


-2.67±0.09 


57 ±22 


(4.1±3.2)e-6 


3552 




-3.03±0.14 


<30 


2.3±0.9)e-6 


7677 


0.01±0.90 


-2.81±0.53 


86 ±47 


(2.2±2.0)e-6 


3569 


-1.47±0.14 




227 ± 114 


2.9±0.9)e-6 


7684 


-0.58±0.20 


-2.10 


646 ± 257 


(9.1±4.5)e-6 


3869 


-1.20 




172 ± 53 


1.3±0.4)e-6 


7750 


-1.20 




76 ± 14 


(2.2±0.5)e-6 


3875 




-3.13±0.69 


<30 


l.l±0.2)e-6 


7769 


-1.50±0.17 




72 ±26 


(1.2±0.5)e-5 


3893 


-0.65±0.04 




153 ±6.3 


2.0±0.1)e-5 


7781 


-0.66±0.77 




94 ±61 


(2.1±1.9)e-6 


4048 


-0.63±0.08 


-2.87±0.45 


323 ± 42 


l.l±0.4)e-5 


7838 


-0.95±0.35 




239 ± 153 


(3.6±2.6)e-6 


4146 


-0.55±0.57 


-3.12±0.51 


85 ±34 


3.2±1.5)e-6 


7845 


-0.09±0.42 




551 ±291 


(4.9±3.5)e-6 


4216 


-1.20 




128 ± 64 


3.5±1.7)e-6 


7989 


-1.52±0.11 


-2.30 


360 ± 216 


(3.7±0.9)e-6 


5417 


-0.71±0.19 


-2.01±0.06 


142 ± 41 


1.5±0.6)e-5 


7998 


-1.26±0.31 




81 ±37 


(3.2±2.0)e-6 



Table A2. Sample of BATSE bursts analyzed in this work. The trigger number and the spectral parameters of the fit of the time integrated spectrum are 
reported. In the last column we report the bolometric fluence obtained by integrating the best fit spectrum. For those bursts whose spectrum allows only to 
set a lower/upper limit on £^°j!|]j we report the BATSE catalogue fluence (i.e. >25 keV). When is not possible to constrain the value of ex we performed the 
spectral fit by fixing a to an appropriate value. These values are reported in table without errors. 
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GRB 


z 


a 


/3 


Peak Flux" 


Range 
keV 


^iso 

erg/s 


-Epcak 

keV 


Ref 


970228 


0.695 


-1.54 [0.08] 


-2.5 [ 0.4 ] 


3.7e-6 [ 0.8e-6 ] 


40-700 


9.1e51 [2.18e51 ] 


195 [ 64 ] 


1 


970508 


0.835 


-1.71 [0.1] 


-2.2 [ 0.25 ] 


7.4e-7 [ 0.7e-7 ] 


50-300 


9.4e51 [ 1.25e51 ] 


145 [ 43 ] 


3 


970828 


0.958 


-0.7 [ 0.08 ] 


-2.1 [0.4] 


5.9e-6 [ 0.3e-6 ] 


30-l.e4 


2.51e52[7.7e51] 


583 [117] 


1 


971214 


3.42 


-0.76 [0.1] 


-2.7 [1.1] 


6.8e-7 [ 0.7e-7 ] 


40-700 


7.21e52[ 1.33e52] 


685 [ 133 ] 


1 


980326 


1.0 


-1.23 [0.21 ] 


-2.48 [0.31] 


2.45e-7[0.15e-7] 


40-700 


3.47e51 [ l.eSl ] 


71 [36] 


4 


980613 


1.096 


-1.43 [0.24] 


-2.7 [ 0.6 ] 


1.6e-7[0.4e-7] 


40-700 


2.e51 [6.7e50] 


194 [ 89 ] 


4 


980703 


0.966 


-1.31 [0.14] 


-2.39 [ 0.26 ] 


1.6e-6[0.2e-6] 


50-300 


2.09e52[4.86e51] 


499 [ 100 ] 


1 


990123 


1.600 


-0.89 [ 0.08 ] 


-2.45 [ 0.97 ] 


1.7e-5 [0.5e-5] 


40-700 


3.53e53 [ 1.23e53 ] 


2031 [ 161 ] 


1 


990506 


1.307 


-1.37 [0.15] 


-2.15 [0.38] 


18.6 [0.1 ] 


50-300 


4.18e52[ 1.33e52] 


653 [ 130 ] 


1 


990510 


1.619 


-1.23 [0.05] 


-2.7 [ 0.4 ] 


2.5e-6 [ 0.2e-6 ] 


40-700 


6.12e52[1.07e52] 


423 [ 42 ] 


1 


990705 


0.843 


-1.05 [0.21] 


-2.2 [0.1] 


3.7e-6[0.1e-6] 


40-700 


1.65e52[2.77e51] 


348 [ 28 ] 


1 


990712 


0.433 


-1.88 [0.07] 


-2.48 [ 0.56 ] 


4.1 [0.3] 


40-700 


7.46e50[ 1.91e50] 


93 [ 15 ] 


1 


991208 


0.706 






1.85e-5 [0.06e-5] 


20-l.e4 


4.32e52 [ 0.38e52 ] 


313[31] 


2 


991216 


1.02 


-1.23 [0.13] 


-2.18 [0.39] 


67.5 [ 0.2 ] 


50-300 


1.13e53[3.75e52] 


642 [ 129 ] 


1 


000131 


4.50 


-0.69 [ 0.08 ] 


-2.07 [ 0.37 ] 


7.89 [ 0.08 ] 


50-300 


1.41e53[5.59e52] 


714 [142] 


1 


000210 


0.846 






2.42e-5[0.15e-5] 


20-l.e4 


8.78e52[l.le52] 


753 [ 26 ] 


2 


000418 


1.12 






2.8e-6 [ 0.4e-6 ] 


20-l.e4 


2.e51 [4.8e50] 


284 [ 21 ] 


2 


000911 


1.06 


-1.11 [0.12] 


-2.32 [ 0.41 ] 


2.0e-5 [ 0.2e-5 ] 


15-8000 


1.65e53[2.89e52] 


1856 [371.] 


1 


000926 


2.07 






1.5e-6[0.26e-6] 


20-l.e4 


4.73e52[1.3e52] 


310. [20. ] 


2 


010222 


1.48 






5.7e-7 [ 0.32e-7 ] 


20-l.e4 


7.87e51 [4.51e50] 


766 [ 30. ] 


2 


010921 


0.45 


-1.6 [0.1] 




9.2e-7[1.4e-7] 


20-l.e4 


7.33e50[ 1.33e50] 


129. [26. ] 


2 


011211 


2.140 


-0.84 [ 0.09 ] 




5.0e-8[l.e-8] 


40-700 


3.17e51 [0.32e51] 


185 [ 25 ] 


1 


020124 


3.198 


-0.87 [0.17] 


-2.6 [ 0.65 ] 


9.4 [1.8] 


2.-400 


5.12e52[2.03e52] 


390 [113] 


1 


020405 


0.695 


-0.0 [ 0.25 ] 


-1.87 [0.23] 


5.e-6 [ 0.2e-6 ] 


15-2000 


1.38e52[7.83e50] 


617 [171 ] 


5 


020813 


1.255 


-1.05 [0.11 ] 




32.3 [2.1] 


2-400 


2.58e52[2.4e51 ] 


478 [ 95 ] 


1 


0208 19B 


0.41 


-0.9 [0.15] 


-2.0 [ 0.35 ] 


7.e-7 [ 0.7e-7 ] 


25-100 


1.49e51 [3.23e50] 


70. [21.] 


7 


020903 


0.25 


-1.0 [0.0] 




2.8 [0.7] 


2-400 


6.7e48 [ 0.26e48 ] 


3.37 [1.79] 


6 


021004 


2.335 


-1.0 [0.2] 




2.7 [0.5] 


2-400 


4.6e51 [0.12e51] 


267 [117] 


6 


021211 


1.01 


-0.85 [ 0.09 ] 


-2.37 [ 0.42 ] 


30 [ 2 ] 


2-400 


7.13e51 [9.9e50] 


94 [ 19 ] 


1 


030226 


1.986 


-0.9 [0.2] 




2.7 [0.6] 


2-400 


8.52e51 [2.23e51] 


290 [ 63 ] 


1 


030328 


1.520 


-1.14 [0.03] 


-2.1 [0.3] 


11.6 [0.9] 


2-400 


l.le52[ 1.55e51 ] 


328 [ 35 ] 


1 


030329 


0.169 


-1.32 [0.02] 


-2.44 [ 0.08 ] 


451 [25] 


2-400 


1.91e51 [2.37e50] 


79 [3] 


1 


030429 


2.656 


-1.1 [0.3] 




3.8 [0.8] 


2-400 


7.6e51 [ 1.47e51 ] 


128 [ 37 ] 


1 


040924 


0.859 


-1.17 [0.05] 




2.6e-6 [ 0.3e-6 ] 


20-500 


6.1e51 [l.leSl] 


102 [ 35. ] 


1 


041006 


0.716 


-1.37 [0.14] 




1.0e-6[0.1e-6] 


25-100 


8.65e51 [ 1.36e51 ] 


108 [ 22 ] 


1 


050126 


1.29 


-0.75 [0.44 ] 




0.698 [ 0.07 ] 


15-150 


1.12e51 [0.25e51] 


263 [110] 


8 


050223 


0.5915 


-1.5 [0.42] 




0.7 [0.1 ] 


15-150 


1.43e50[0.2e50] 


110 [54] 


8 


050318 


1.44 


-1.34 [0.32] 




3.2 [0.3] 


15-150 


5.11e51 [0.8e51 ] 


115 [27] 


8 


050401 


2.9 


-1.0 [0.0] 


-2.45 [0.0 ] 


2.45e-6[0.12e-6] 


20-2000 


2.03e53[0.1e53] 


501 [117] 


9 


0504 16A 


0.653 


-1.01 [0.0] 


-3.4 [0.0 ] 


5.0 [0.5] 


15-150 


9.3e50 [ 0.9e50 ] 


28.6 [ 8.3 ] 


10 


050505 


4.27 


-0.95 [0.31 ] 




2.2 [ 0.3 ] 


15-350 


5.65e52 [ 0.8e52 ] 


661. [245] 


11 


050525A 


0.606 


-0.01 [0.11] 




47.7 [1.2] 


15-350 


9.53e51 [2.5e51 ] 


127 [5.5] 


12 


050603 


2.821 


-0.79 [0.06 ] 


-2.15 [0.09] 


3.2e-5 [ 0.32e-5 ] 


20-3000 


2.13e54[0.22e54] 


1333 [ 107 ] 


13 


050803 


0.422 


-0.99 [0.37 ] 




1.5 [0.2] 


15-350 


1.31e50[2.6e49] 


138 [48] 


14 


050814 


5.3 


-0.58 [0.56 ] 




1.0 [0.3] 


15-350 


3.0e52[5.6e51] 


339 [ 47 ] 


15 


050820A 


2.612 


-1.12 [0.14] 




1.3e-6[0.13e-6] 


20-1000 


9.1e52[6.8e51] 


1325 [ 277 ] 


16 


050904 


6.29 


-1.11 [0.06] 


-2.2 [0.4 ] 


0.8 [ 0.2 ] 


15-150 


I.le53[3.9e52] 


3178 [ 1094.] 


17 


050908 


3.344 


-1.26 [0.48] 




0.7 [0.1] 


15-150 


8.29e51 [ 1.3e51 ] 


195 [ 36 ] 


18 


050922C 


2.198 


-0.83 [0.26 ] 




4.5e-6 [ 0.7e-6 ] 


20-2000 


1.9e53 [2.3e51] 


417 [118] 


19 


051022 


0.80 


-1.176 [0.038] 




l.e-5[0.13e-5] 


20-2000 


3.57e52[2.7e51 ] 


918 [63] 


20 


051109A 


2.346 


-1.25 [0.5] 




5.8e-7 [ 2.e-7 ] 


20-500 


3.87e52[3.8e51 ] 


539 [381] 


21 


060115 


3.53 


-1.13 [0.32] 




0.9 [0.1] 


15-150 


1.24e52[2.0e51] 


288 [ 47 ] 


22 


060124 


2.297 


-1.48 [0.02] 




2.7e-6 [ 0.8e-6 ] 


20-2000 


1.42e53[1.35e51] 


636 [ 162 ] 


23 


060206 


4.048 


-1.06 [0.34] 




2.8 [0.2] 


15-150 


5.57e52[9.0e51 ] 


381 [98] 


24 


060210 


3.91 


-1.12 [0.26] 




2.7 [0.3] 


15-150 


5.95e52[8.0e51] 


575 [186] 


25 


060218 


0.0331 


-1.622 [0.16] 




l.e-8[0.1e-8] 


15-150 


1.34e47[0.3e47] 


4.9 [0.3] 


26 


060223A 


4.41 


-1.16 [0.35] 




1.4 [0.2] 


15-150 


3.27e52[5.5e51 ] 


339 [ 63 ] 


27 


060418 


1.489 


-1.5 [0.15] 




6.7 [0.4] 


15-150 


1.89e52 [ 1.59e51 ] 


572 [114] 


28 


06051 OB 


4.9 


-1.53 [0.19] 




0.6 [0.1] 


15-150 


2.26e52 [ 1.78e51 ] 


575 [ 227 ] 


29 



Table A3, continue.. 
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GRB 


^ 


a 


/3 


Peak Flux" 


Range 
keV 


erg/s 


-Epcak 

keV 


Ref 


060522 


5.11 


-0.7 [0.44 ] 




0.6 [0.1] 


15-150 


2.0e53[3.7e51] 


427 [ 79 ] 


30 


060526 


3.21 


-1.1 [0.4] 


-2.2 [0.4 ] 


1.7 [0.1] 


15-150 


1.72e52[3.1e51] 


105.2[21.1] 


31 


060605 


3.78 


-1.0 [0.44] 




0.5 [0.1] 


15-150 


9.5e51 [1.5e51] 


490 [251] 


32 


060607A 


3.082 


-1.09 [0.19] 




1.4 [0.1] 


15-150 


2.0e52[2.7e51] 


575 [ 200 ] 


33 


060614 


0.125 






11.6 [0.7] 


15-150 


5.3e49[ 1.4e49] 


55 [ 45 ] 


34 


060707 


3.43 


-0.73 [0.4 ] 




1.1 [0.2] 


15-150 


1.4e52[2.8e51] 


302 [ 42 ] 


35 


060714 


2.711 


-1.77 [0.24] 




1.4 [0.1] 


15-150 


1.42e52[l.e51] 


234 [ 109 ] 


36 


060814 


0.84 


-1.43 [0.16] 




2.13e-6[0.35e-6] 


20-1000 


I.e52[l.e51 ] 


473 [ 155 ] 


37 


060904B 


0.703 


-1.07 [0.37] 




2.5 [0.1] 


15-150 


7.38e50[1.4e50] 


135 [41] 


38 


060906 


3.686 


-1.6 [0.31] 




2.0 [ 0.3 ] 


15-150 


3.55e52[3.9e51] 


209 [ 43 ] 


39 


060908 


2.43 


-0.9 [0.17] 




3.2 [ 0.2 ] 


15-150 


2.6e52[4.6e51] 


479 [110] 


40 


060927 


5.6 


-0.93 [0.38 ] 




2.8 [0.2] 


15-150 


1.14e53 [2.0e52] 


473 [116] 


41 


061007 


1.261 


-0.7 [0.04 ] 


-2.61 [0.21 ] 


1.95e-5 [0.28e-5] 


20-le4 


1.74e53 [2.45e52] 


902 [ 43 ] 


42 


061121 


1.314 


-1.32 [0.05 ] 




1.28e-5[0.17e-5] 


20-5000 


1.41e53[1.5e51] 


1289 [ 153 ] 


43 


061126 


1.1588 


-1.06 [0.07] 




9.8 [0.4] 


15-150 


3.54e52[3.0e51] 


1337 [410] 


44 


061222B 


3.355 


-1.3 [0.37] 




1.5 [0.4] 


15-150 


1.82e52[2.75e51] 


200 [ 28 ] 


45 


070125 


1.547 


-1.1 [0.1] 


-2.08 [0.13] 


2.25e-5 [ 0.35e-5 ] 


20-l.e4 


3.24e53 [ 5.e52 ] 


934 [ 148 ] 


46 


070508 


0.82 


-0.81 [0.07] 




8.3e-6[l.le-6] 


20-1000 


3.3e52[3.9e51] 


342 [ 15 ] 


47 


071003 


1.100 


-0.97 [0.07 ] 




1.22e-5[0.2e-5] 


20-4000 


8.4e52 [ 1.5e51 ] 


1678 [231] 


48 


071010B 


0.947 


-1.25 [0.6] 


-2.65 [0.35 ] 


8.92e-7 [ 3.7e-7 ] 


20-1000 


6.4e51 [5.3e49] 


101 [23] 


49 


071020 


2.145 


-0.65 [0.3 ] 




6.04e-6[2.1e-6] 


20-2000 


2.2e53[9.6e51] 


1013 [ 205 ] 


50 


071117 


1.331 


-1.53 [0.15] 




6.66e-6[1.8e-6] 


20-1000 


I.e53[7.e51] 


648 [318] 


51 


0803 19B 


0.937 


-0.82 [0.01 ] 


-3.87 [0.8 ] 


2.17e-5[0.21e-5] 


20-7000 


9.6e52[2.3e51] 


1261 [25] 


52 


0803 19C 


1.95 


-1.2 [0.1 ] 




3.35e-6 [ 0.74e-6 ] 


20-4000 


9.5e52[ 1.2e51] 


1752 [ 505 ] 


53 



Table A3. "Peak Fluxes are in erg/s/cm^ or photons/s/cm^ . Reference for the Peak Flux (or Luminosity): 1) Firmani et al. 2006; 2) Ulanov et al. 2005 (Liso 
computed as [P/F] [1 + zji^iso); 3) Jimenez et al. 2001 ; 4) Amati et al. 2002; 5) Price et al. 2003; 6) Sakamoto et al. 2005; 7) Hurley et al. GCN 1507; 8) GRB 
BAT on line table (http://swift.gsfc.nasa.gOv/docs/swift/archive/grb_table/l; 9) Golenetskii et al. 2005f, GCN 3179; 10) Sakamoto et al. 2006a 11) Hullinger et 
al. 2005; 12) Blustin et al. 2006; 13) Golenetskii et al. 2005g; 14) Parson et al. 2005; 15) Tueller et al. 2005; 16) Cenko et al., 2006; 17) Sakamoto et al. 2005b; 
18) Sato et al. 2005; 19) Golenetskii et al. 2005h; 20) Golenetskii et al. 20051; 21) Golenetskii et al. 2005J; 22) Barbier et al. 2006a; 23) Golenetskii et al. 
2006J; 24) Palmer et al. 2006a; 25) Sakamoto et al. 2006b; 26) Campana et al., 2006; 27) Cummings et al. 2006a; 28) Cummings et al. 2006b; 29) Barthelmy 
et al. 2006; 30) Krimm et al. 2006a; 31) Markwardt et al. 2006a; 32) Sato et al. 2006a; 33) Tueller et al. 2006; 34) Mangano et al. 2007; 35) Stamatikos et 
al. 2006a; 36) Krimm et al. 2006b; 37) Golenetskii et al. 2006k; 38) Mai-kwardt et al. 2006b; 39) Sato et al. 2006b; 40) Palmer et al. 2006b; 41) Stamatikos 
et al. 2006b; 42) Golenetskii et al. 20061; 43) Golenetskii et al. 2006m; 44) Krimm et al. 2006c; 45) Barbier et al. 2006b; 46) Golenetskii et al. 20071; 47) 
Golenetskii et al. 2007m; 48) Golenetskii et al. 2007n; 49) Golenetskii et al. 2007o; 50) Golenetskii et al. 2007p; 51) Golenetskii et al. 2007q; 52) Golenetskii 
et al. 2007r; 53) Golenetskii et al. 2007s. 
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